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A METHOD FOR SOLVING THE NONSIMILAR LAMINAR 


BOUNDARY- LAYER EQUATIONS INCLUDING 
FOREIGN GAS INJECTION 

By Joseph G. Marvin and Yvonne S. Sheaffer 
Ames Research. Center 


SUMMARY 


A numerical method for solving the nonsimilar boundary- layer equations, 
including binary gas injection, is developed for nonreacting gases. The con- 
tinuity, momentum, energy, and species-concentration equations and boundary 
conditions are reduced to a set of linear algebraic equations in terms of the 
dependent variables only by appropriate application of implicit finite differ- 
ence expressions. To solve the equations, profiles of the unknowns at the 
initial station are specified, together with certain boundary conditions along 
the body surface and boundary- layer edge. The initial profiles are obtained 
through solution, by the same method, employing the similarity form of the con- 
servation equations that apply exactly at the starting points for the problems 
considered. 

Comparisons of the present solutions with exact similarity solutions for 
air and helium injection are presented to demonstrate the accuracy of the 
method for evaluating the dependent variables along the initial data line. 
Comparisons of the present solutions with nonsimilar solutions obtained from 
other numerical methods and comparisons with experimental data for air, argon, 
and helium injection are presented to evaluate the accuracy and advantages of 
the method. 


INTRODUCTION 


At the inception of this investigation, wind-tunnel tests were planned 
that involved foreign gas injection through a porous model. To meet the 
requirements for data correlations, it was necessary to calculate various 
boundary- layer parameters that required solutions to the nonsimilar boundary- 
layer equations including foreign gas injection. A computer program utilizing 
the Smith-Clutter approach was available for solving the nonsimilar equations 
(see, e.g., refs. 1 and 2). While experience has shown this program was use- 
ful for studying the effects of pressure gradient without injection (see 
refs. 3 and 4), computer time and convergence limitations made it impractical 
for the present application. Therefore, a finite difference method for 
solving the nonsimilar equations including binary gas injection was developed 
to circumvent these limitations. 



Finite difference schemes can be categorized as explicit or implicit. 

The explicit schemes are conditionally stable; for most applications, the Ax 
step size becomes impractically small and resulting computing times exces- 
sively long (see ref. 5). The implicit schemes, on the other hand, are 
usually unconditionally stable. Fliigge-Lotz and Blottner (ref. 6) took advan- 
tage of this fact to demonstrate the utility and speed of an implicit scheme 
for solving the boundary- layer equations. Blottner later applied this tech- 
nique (ref. 7) to solve the complicated nonequilibrium boundary- layer problem. 
Although the method described in reference 6 has many advantages, it was not 
applied extensively until recently (e.g., refs. 8 and 9). 

The present work uses the basic implicit scheme described in reference 6, 
but here the basic equations are transformed to new coordinates with finite 
boundaries by means of a transformation described in reference 10. The asymp- 
totic nature of the transformation allows the choice of a constant Ay step 
size within finite limits (say, between 0 and 1), thus eliminating the neces- 
sity of searching for the effective boundary- layer edge through additional 
iteration. Subsequent to the coordinate transformation, the equations are put 
in finite difference form so that the unknowns appear linearly. The resulting 
set of linear algebraic equations and a stream- function equation are solved on 
a digital computer after specifying streamwise boundary conditions along the 
body surface and boundary- layer edge and profiles of the unknowns at the 
initial streamwise location. The initial profiles are obtained by the same 
method employing the similarity form of the equations that apply exactly at 
the streamwise starting locations for the problems considered. 

This report presents the solution method and demonstrates, by example, 
its capabilities and accuracy. Although the original application of the 
method was specifically directed toward solutions that included binary 
injection, the method is applicable to many other boundary- layer problems. 


SYMBOLS 


A n ,B n ,C n coefficient matrices given in equation (23) 

C mass concentration 

Cp frozen specific heat at constant pressure, see equation (18) 

specific heat at constant volume 
D n solution vector (column matrix) in equation (23) 

D* total displacement thickness with mass addition, see equation (39) 

D 12 binary diffusion coefficient 
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defined by equations (28) and (29) 

defined by equation (4) 

„ f n u , 

stream function, I — dn 
J U e 

dummy variable 

edge boundary-condition vector 

boundary-condition coefficient matrix 

wall boundary-condition vector denoting column matrix 

vector denoting column matrix given by equation (28) 

weighting factor in transport coefficients, see equation (22) 

thermal conductivity 

length of nonporous tip 

dummy variable 

molecular weight 

mesh point in x direction 

denotes number of mesh points in y direction 

mesh point in y direction, n = 1,2, . . . , N 

coefficient in polynomial expressions for properties, see 
equation (15) 

pressure 

Prandtl number, ~ 

9 P k 


heat-transfer rate 
dummy variable 
universal gas constant 
normal distance from axis of 
dimensionless distance along 


symmetry to body surface 
surface in flow direction -r- 

H 


Schmidt number. 
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normalized temperature, ~- 

r o 


temperature 

velocity in streamwise direction 
velocity normal to streamwise direction 
mole fraction 

streamwise distance from tip of body 

normal distance from body surface 

transformed normal distance defined by equation (7) 

stretched normal distance, a(l - y) 

scale factor in equation (7) 


, . 2£ du e 

pressure gradient parameter, 

displacement thickness given by equation (38) 

wall boundary conditions on species (eqs. (26) and (30)) 

transformed normal coordinate (eq. (5)) 

momentum thickness given by equation (40) 

viscosity 

transformed streamwise coordinate (eq. (6)) 
density 

vector denoting column matrix for unknowns at (m + 1), see 
equation (23) 


shear stress defined by equation (33) 
dimensionless density-viscosity product, 


py 

Pe^e 


step size in y direction 



Subscripts 


a,b,c,d,e,f mesh points for finite difference scheme (see sketch (b)) 
e boundary- layer-edge value 

i species i; i = 1 is free stream, i = 2 is injectant 

I initial value 

k coefficient number in polynomial for thermodynamic and transport 

properties 

V evaluated at arbitrary length, x = Z 

0 stagnation-point value 

NS nonsimilar value 

S similar value 

y,Ti,s,5 partial differentiation with respect to the subscript 

W wall value 

1 free-stream species 

2 injected species 

Superscripts 

v iteration number 

n geometric factor; n = 0 is two-dimensional and n = 1 is 

axisymmetric 

k polynomial degree 


ANALYSIS 


The nonsimilar, laminar boundary- layer equations are given in this 
section, along with the equations for evaluating the thermodynamic and trans- 
port properties of the binary gas systems considered. The method of solution, 
involving a transformation and a linearization of the boundary- layer equations, 
is presented and equations are given for determining shear stress, heat trans- 
fer, and other important boundary- layer parameters. 


5 



' ■ l III HI 


Boundary- Layer Equations 

The mass, momentum, energy, and species conservation equations for a 
binary gas system in the familiar Levy-Lees variables (e.g., ref. 2) are: 


Momentum 


(♦F n l n + £F n + B(P e / P - F2) = 25(FF 5 - P n f 5 ) 


Energy 

1 l ♦ r t\ ♦ £T * _i_ Ic .C \c T 

C^^T7 C P T nj n fT n ScCp ( C Pj P 2 ) \ " c p t 0 n c p t 0 p 


Species (binary system) 


'± Cl ) * fC ln - 2 5 (fC 1{ - C, n f 5 ) 


where 


f = ± 
n u e 

rI1 Pe u e cy _e_ d y 

m Jo Pe 


5=1 P„u-V*r 2n dx 


The boundary conditions are 
n = 0 


/2f (cA (pv) Sc 
\ / w w v 


Pw% u e r 


f (0) = (pv) w r n dx 

m Jx i 
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n •> 


oo 


F 1 
T ^ T e 


The boundary conditions at n = 0 are based on the assumption that the net 
mass flow of the free-stream gas at the wall is zero (e.g., ref. 2). 


Transformation to Finite Coordinates 

It is convenient from a computational viewpoint to transform the above 
equations to a new system of coordinates wherein the indefinite limit of inte- 
gration on n is replaced by a definite limit (e.g., ref. 10). Let 

y = 1 - e" a F (7) 

For a system of finite-difference equations with a fixed number of y nodal 
points and a fixed interval Ay, a is used as a scaling factor to provide an 
optimum distribution of nodal points across the boundary layer. A certain 
amount of experience is required to achieve the optimum results, as will be 
discussed subsequently. 

The following transformation identities result 

t 1 ., = 0(1 - y)( ) y 

( - “ 2 (1 - yH(i - y)C 1 yy - ( )y] (8) 

Employing these transformation identities, equations (1) through (4) 
become 


Momentum 


a 2 (l - y) 2 <)>yFy + <t>a 2 (1 - y) 2 Fyy “ 4>a 2 (l - y)F y + a(l - y)fF y + £(p e /p - F 2 ) 


= 2C[FF 5 - a(l - y)Fy£ 5 ] 


Energy 


^a 2 (l- y)2(^E) yT y + ^a 2 (l _ y) [Cl - y)T yy - T y ] + af(l - y)T y 

+ ( c p, " c P 2 ) a2(1 - y)2c i y T y " c4; ~t F + a2(1 ' y)2? y 2 

= 2^[FT ? - a(l - y)T y f c ] 


O) 


(10) 
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Species 


“ 2( 1 " y)2C l y(jk)y + SF “ 2(1 - y)2 % y - SF ° 2(1 - y)C ly + a£(1 - y)C ly 

= 2?[fC 1 - a (1 - y)C^fJ (11) 


where 


J y p 

-7— v dy + f (0) 

0 a(l - y) 3 

The boundary conditions become for y = 0 

-\ 

F = 0 
T = T w 

(pv) Sc 

w W l 

ap w^w u e r11 

f (0) = — f X (pv) x n dx 

m J w J 



and for y = 1 


(12) 


(13) 


F = 
T = 


C 


1 ~ 



(14) 


Thermodynamic and Transport Properties 

To complete the system of equations, the thermodynamic and transport 
properties are expressed as functions of the dependent variables C 1 and T in 
the manner described in reference 2. The properties of the individual gas con- 
stituents were expressed as polynomials in terms of temperature; for example. 



(15) 


A similar expression was used for and Di2- Values of the coefficients, P^, 
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used to obtain y^, and D — for the gases considered are given in 

reference 2, along with a brief description of the assumptions regarding 
intermolecular potentials. 

The thermal conductivity of each constituent was obtained from Eucken*s 
relation 



where 



(16) 


(17) 


Mixture values of the thermodynamic and transport properties were 
obtained from the equations : 

C P = S 2 + c i( c p 1 ' c p 2 ) 

p r M i M i ] 

p RT [ M[ + C J 


(18) 


(19) 


Mixture values of the transport properties were obtained from the approximate 
combining equations taken from reference 11: 


k 


1 + 1.065K 12 ^ 
A l 


1 + 1.065K 21 


( 20 ) 


y 




l + k 12 


X2 

Xi 


+ 


y 2 


i + k 21 


Xl 

x 2 


where 


k 12 = 


[1 + (y 1 /y 2 ) 1/2 (^ 2 /M 1 ) I/4 ] 2 
2 3/2 (l + M 1 /M z y /Z 


( 21 ) 


( 22 ) 


and where K 2 i is obtained from equation (22) with subscripts 1 and 2 
interchanged. 
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Method of Solution 


To obtain a numerical solution, equations (9) through (11) are 
approximated by a set of implicit finite difference equations of the Crank- 
Nicholson type and the ensuing tridiagonal matrix equations A53 = D, solved on 
a computer by the algorithm given in reference 12. To use the solution algo- 
rithm, the equations must be linearized, values for the unknowns at the 
initial streamwise location must be specified along with streamwise boundary 
conditions at the wall and the boundary- layer edge. In the subsequent por- 
tions of this section, the difference equations, the boundary conditions, and 
the solution are discussed. 

The mesh-point diagram for the Crank-Nicholson scheme is shown below. 

The finite-difference equations used to assure linear equations when marching 

from line (m) to (m + 1) are given in appendix A. 

Substituting the difference equations into the 
momentum, energy, and species equations (eqs. (9) - 
(11)) results in 3(N - 2) linear equations. In 
matrix equation form, one obtains 

An%- 1 + ®n w n + Qn (ii n+ i = n = 2,3, . . . , N - 1 

(23) 

where, for example, the vectors and coefficient 
matrices have the following form 



m m + l 

Sketch (a) 



or n,m+i/2 or n,m+l/2 

As shown in appendix A, the matrices B n and C n are similar in form to A n , 
except B 2 3 = 0. The matrices A n , B n , and C n , and the solution vector S n , 

are considered known and are evaluated at (n, m) or (n, m + 1/2) depending on 
the iteration described subsequently. Equations for the matrix elements and 

vector S n are given in appendix A. 


To complete the system of equations, values of the dependent variables 
u> n at the two boundaries must be specified. To be consistent with the 
solution algorithm in reference 12, these were written as 


y = 0 = Hu) 2 + H 

y = 1 = I 


(24) 
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Special attention was given to the boundary condition on wall concentration to 
arrive at a form consistent with equation (24) . The condition was obtained 
by expressing the concentration at n = 2 in finite difference form as a 
Taylor series about n = 1 and retaining three terms. 

(Ci) 2 = CjCAy) = CjCO) + Ay [Ci(0)] y + [C x C0)] yy (25) 


The terms [Ci(0)] y and [Ci(0)] yy are obtained from the boundary condition, 

equation (13), and equation (11) evaluated at y = 0. After some manipulation, 
equation (25) is rewritten as 


where 


Ci(0) 


Cl (Ay) 

1 + e NS 


(26) 


e NS 


Sit Sc w (pv) w 

— — Ay + 

a(pvO w u e r n 




/2£ (pv) Sc 
v ^w w 

a(pvi) w u e r n 


Using equation (26) together with equations (13) and (14) , the elements of the 
matrix coefficient and vectors in equations (24) become 


H = 


0 

0 


0 

0 

-1 


1 + e NS 


h = 


W 


g = 


n,m+i 


u 


71,111 

or n,m+i/2 


Equations (23), together with the boundary conditions, equations (24), 
form a set of linear equations which can be solved on a digital computer once 
profiles of the dependent variables at an initial streamwise location are 
specified. The algorithm for solution, taken from reference 12, is 


n 


= -E^w. 


n^n+i 


+ J 


n 


1 < n < N - 1 


(27) 


where 

E„ ' On - VW-'C,, 

3 n - < B n - VW* 2 s n < N - 1 

and where 

Ei = H E n = 0 

Jl = h Jn = g 


(28) 


(29) 
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The computation proceeds with evaluation of equations (28) from the wall 
to the boundary- layer edge using matrix elements based on values of the depen- 
dent variables and appropriate derivatives evaluated at m. Values of the 
dependent variable 5 n at (m + 1) are determined in reverse order using equa- 
tions (27) and (29). New values of the stream function f are obtained by 
integrating equation (12) . Values of the streamwise derivatives of f are 
redetermined and the new values for f and u n are averaged over the inter- 
val with corresponding values at m to form new matrix elements at 
(n, m + 1/2). The solution was iterated until 


|[F(0)]j; +l - [H0)]y 


[HO)] 

(For all computations made thus far, this criterion was sufficient to assure 
convergence of temperature and concentration profiles.) 


v+i 


< 0.0005 


To start the solutions, values of the dependent variables 5 n and the 
stream function f at the initial station ?j, must be specified. From an 
appropriate starting point these values can be obtained by solving the simi- 
larity form of the boundary- layer equations obtained by setting the right-hand 
sides of equations (9) through (11) to zero (i.e., visualize a stagnation 
point, or a point along the surface of a cone with uniform free-stream condi- 
tions and f(0) = 0 or f(0) ~ l//x). The similarity form of the equations 
were solved by a method analogous to that described above. The mesh-point 
diagram and finite difference approximations are given in appendix B. Substi- 
tuting the difference equations from appendix B into equations (9) through 
(11) with the right-hand sides set to zero also results in 3(N - 2) linear 
equations of the form given by equation (23) . The equations for the coeffi- 
cient matrix elements and vector B n are given in appendix B. The boundary 
conditions and solution algorithm remain unchanged, except that the term 
e^jg in the concentration boundary condition becomes 


/2f Sc (pv) 
w VK w 

[ Ay * ^2“ 

r sc 
1 w 

/<p \ 

w 

Sc 

e s - 

a(pu) u r n 

1 <p ' 
w 

1 Sc J 
\ w/. 

<f> i 

w 

Y 


and the boundary condition on the stream function becomes 



(30) 


f(0) 


(pv) w 

(pv) e u e rn 


(31) 


The similarity solution is started by assuming linear profiles of the depen- 
dent variables across the boundary layer to evaluate the matrix elements in 
equations (28). New values of the dependent variables are obtained from equa- 
tions (27) and (29) and integration of equation (12) . The new variables are 
then averaged with corresponding values of the variables from the previous 
iteration to form new matrix elements for equations (28). The solution was 
iterated until 


[F (0) ]y +1 - [F(O)]^ 

[F(o)]j ; +1 


0.0005 
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Boundary- Lay er Parameters 


In addition to using the solutions of the foregoing boundary -layer 
equations to obtain profiles of F, T, and C at each station along a body, 
other quantities can be computed as follows: 

Skin friction .- The defining equation for wall skin friction is 

T w = (32) 

Transforming equation (32) to the present coordinate system gives 

x = a ^e 2 PwP W r n /p \ (33) 

/2? V /w 

The derivative (F y ) was obtained by means of a three-point derivative 
v 7/ w 

equation and the converged solution for F 

( F y) w = tV(" 3F » + 4F **«y ■ F »* 2 ay) (54) 

Heat transfer.- The heat transfer to the wall was obtain from 


-tw 


= -to^w(T^) w 


atoPwkw rn u e 

/2l 


w 


(35) 


where (Ty) w was obtained from the converged solution for T in an equation 

similar in form to equation (34). At a stagnation point, equation (35) 
reduces to 

J (du /dx) Q 

2n (T y } w (36) 

r e e 

Displacement and momentum thicknesses .- Certain useful integrals of the 
profiles across the boundary layer can also be evaluated. The usual displace- 
ment thickness is given by 


6 * 



pu 

p e u e 



(37) 


In the present coordinate system this equation becomes 


6* 


/2T f ye (Pe/P)~ P . 

Pe u e rn Jo a(1 " y) 7 


(38) 
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where y e is the boundary- layer height at which F = 0.995. It should be 
noted that the thickness given by 6* is not the displacement thickness 
required to describe physically the fvee-stveam mass flow entrained in the 
boundary layer with injection. To obtain that thickness, as shown in refer- 
ence 13, the appropriate expression is 


D* = <5* 


if. 

r n J 


(pv) 


p e u e 


& r n dx 


(39) 


The momentum thickness is given by the equation 
e - f e (l - — ) dy = S— f } 

1 PeUe ^ Ue ' w n J 0 


F Cl - F) 
a(l - y) 


dy 


(40) 


RESULTS AND DISCUSSION 




(a) Velocity profile. 

(b) Temperature profile. 

Figure 1.- Comparison of the present solution 
with the solution of Low, reference 14. 


The foregoing equations were 
programmed in Fortran IV for solu- 
tion on an IBM 7094 computer. The 
method of solution was verified by 
comparisons with numerical solu- 
tions obtained from other sources 
and by comparisons with experi - 
mental data. Some examples of 
these comparisons are given below 
to illustrate the capabilities and 
advantages of the method. 

Similarity Solutions 

As noted in the analysis 
section, a special finite- 
difference formulation was used to 
solve the similarity form of the 
boundary- layer equations. To ver- 
ify this approach, a comparison 
with exact similarity solutions 
was made. Typical examples are 
discussed below. 

Figure 1 shows a comparison 
with the velocity and temperature 
profiles from reference 14. Pro- 
files, with and without air injec- 
tion, are shown. (The linear 
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viscosity law reported in reference 14 was also employed in the present 
solutions to eliminate transport property differences.) Results from the pres- 
ent method were obtained for several values of the scale factor a using 21 
nodal points (N) across the boundary layer. As mentioned previously, the 
choice of a determines the spacing of the nodal points across the layer and 
the choice has some influence on the solution in the outer portion of the 
boundary layer. See, for example, the solution for f(0) = 0 in figure 1. 
Experience has shown values of a < 0.5 are required when injection is being 
considered. However, even for a poor choice of a, the profiles and deriva- 
tives of velocity and temperature near the wall are accurate enough for most 
practical applications. No attempt was made to determine the optimum combina- 
tions of a and N for this or any of the following examples. 

A comparison of the present 
solution with the similarity solu- 
tion from reference 15 for for- 
eign gas injection is shown in 
figure 2. The thermodynamic and 
transport properties given in the 
analysis section were used for 
the present solutions, and these 
were not exactly the same as 
those reported in reference 15. 
Velocity, helium-concentration, 
and temperature profiles agree 
very well, except for differences 
in the concentration near the 
wall. The differences were attri- 
buted to the different transport 
properties employed in the two 
solutions. 


Nonsimilar Solutions 

The present method was 
verified primarily through compar- 
isons with solutions using a modi- 
fied form of the Smith-Clutter 
computer program (see refs. 1 and 
3). In the Smith-Clutter method, 
the partial differential equa- 
tions are reduced to ordinary dif- 
ferential equations in n by 
finite differencing derivatives 
in the £ direction. The equa- 
tions are solved by the so-called "shooting" technique - that is, as a two- 
point boundary-value problem prescribing initial conditions for derivatives 
of velocity and temperature at n = 0, and iteratively changing these condi- 
tions until the edge conditions are satisfied as n -* 00 . 




(a) Velocity and concentration. 

(b) Temperature. 


Figure 2.- Comparison of the present solution 
with the solution of Baron, reference 15. 
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(a) Heating. 

(b) Shear stress. 


Figure 3.- Comparison of nonsimilar solutions for a sphere. 


Figure 3 presents an example comparison of predicted heating rates and 
wall shear along the surface of a hemisphere. The same values of surface 
pressure and boundary- layer- edge velocity and temperature were used as input 
to both computer programs. Agreement between the two methods is excellent, 
as was the agreement between velocity and temperature profiles at various loca 
tions along the hemisphere. As described in the analysis section, initial pro 
files for the present method were obtained by solving the finite difference 
form of equations (9) through (12) at the stagnation point where £ = 0. Com- 
putations were then made for the full set of equations at 92 locations along 
the hemisphere. The Ax step size was constant at a value of 0.005 ft. No 
attempt was made to determine the maximum Ax step size to achieve minimum 
computing time because such optimization would depend on body geometry and 
corresponding boundary- layer- edge conditions and the results would not be gen- 
erally applicable to other examples. However, it is interesting to note that 
the present solutions required about 0.04 minute per body station versus about 
0.7 minute per station for the method of reference 1. 

A more stringent example is given in figure 4. Solutions were obtained 
for a cone with a sharp solid tip followed by a porous surface with uniform 
air injection (see ref. 16). The decrease in skin friction and accompanying 
increase in displacement thickness along the cone surface is compared with 
corresponding values obtained from the computer programs described in refer- 
ences 1 and 16. Boundary- layer- edge conditions and cone pressure were assumed 
constant in all programs. The starting profiles for the present solutions 
were obtained from the similarity form of equations (9) through (12), evalu- 
ated at x = 0.158 ft and for f(0) = 0. Solutions beyond this location were 
obtained from the complete equations with the prescribed value of (pv) w . The 
Ax step size was 0.00025 ft for the first 21 computing stations and then was 
increased to 0.015 ft for the remaining stations. The small initial step size 
was chosen to avoid numerical instabilities due to the step change in (pv) w . 
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For the lowest injection rate, agreement with the other computed solutions is 
very good. For the highest injection rate considered (pv) w /( pu)^ = 2. 37x10” 3 , 

solutions using the "shooting" technique (ref. 1) did not satisfy the deriva- 
tive boundary conditions as n + 00 for x > 0.2 ft; in addition, the program 
would not continue beyond x = 0.325 ft, whereas there was no difficulty in 
proceeding further with the present method and that of reference 16. The dif- 
ficulty encountered by the "shooting" technique is typical of those encount- 
ered by this method when gradients of the dependent variables near the wall 
become small. The distinct advantage of the present method is that values of 
the dependent variables, rather than their slopes, are employed to obtain solu- 
tions. As a result, no difficulty has been encountered with injection prob- 
lems where the skin friction or heat transfer become very small. This 
advantage would also be expected when solving problems with adverse pressure 
gradient approaching incipient separation. It is not known whether the method 
described in reference 16 also has these advantages. In the example above, 
the solution from reference 16 failed to converge at x = 0.65 ft (the next 
computing station), probably because the Ax step size was too large. 

One of the purposes for developing the present finite-difference method 
of solution including foreign gas injection was to assist in data correlations 
of wind-tunnel experiments. A configuration currently under study is similar 
to that used in the previous example; that is, a 5° cone with a solid sharp 
tip followed by a uniformly porous conical surface. The predicted effects of 
gas injection for this test model are given below and compared with some exper- 
imental heat-transfer data taken at Ames Research Center in the 3.5-ft hyper- 
sonic wind tunnel (see ref. 17). 

Figure 5 presents the predicted laminar boundary- layer heat transfer 
along the porous cone surface for three injectant gases and for several inject- 
ant rates. The predictions were obtained assuming uniform boundary- layer-edge 
conditions of Mach number, temperature, and pressure, and using the measured 
injection rate distribution from reference 17. The symbols represent measure- 
ments of the laminar boundary- layer heating rate obtained for each of the 
injectant rates. Only laminar boundary- layer heating data are shown in fig- 
ure 5. Transition occurred and the heating rates increased downstream of the 
locations where these data are plotted. For each injectant gas, increasing 
the injectant rate decreased the heat transfer to a negligible value. For the 
same injectant rate the amount of heating reduction depended on the injectant 
gas molecular weight; compare, for example, the reductions for each injectant 
at (pv^/PqU^ = l.lxltT 3 . The agreement between theory and experiment is good, 
considering the experimental uncertainties in local injection rates and the 
fact that the theory did not account for the changes in boundary- layer- edge 
Mach number, temperature, and pressure resulting from displacement effects due 
to injection. * The displacement effects are particularly important for the 
helium injection tests. 
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L = 0.3125 ft 


f 3 = 0 

... 

Present solution, N = 2I, a = 0.5 

u e = 4029 ft/sec 

SI 


t c = 1460° R 

n 

O l 

Data, ref 17 

t w = 574° R 

A 1 


f e r 144° R 



p e = 26.1 lb/ft 2 



{/?u)co = OJ92 slugs/ft 2 sec 





(a) Air injection. 

(b) Argon injection. 

(c) Helium injection. 

Figure 5.- Comparison of present solution with heat-transfer data on a sharp 5° porous 

cone, reference 17. 


CONCLUDING REMARKS 


An accurate numerical method was developed for solving the nonsimilar 
boundary- layer equations, including binary gas injection. The conservation 
equations in the familiar Levy-Lees coordinates were transformed to new coordi- 
nates having a finite domain. By the appropriate choice of implicit finite 
difference expressions, the governing partial differential equations and 
boundary conditions were reduced to a set of linear algebraic equations in 
terms of the dependent variables only and a single stream function equation. 
Initial values of the dependent variables required to solve the nonsimilar 
equations were determined by applying the same method to solve the similarity 
form of the boundary- layer equations. 

The capabilities and accuracy of the method were demonstrated by 
comparison with other numerical solutions and with experimental data. The com- 
parisons showed that the present method could be accurately applied to a vari- 
ety of examples. The method is ideally suited to solving flow problems where 
skin friction and heat transfer become small, as in the case of boundary- layer 
injection. 
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Additional information regarding the computer program details can be 
obtained by contacting the authors . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, July 16, 1969 
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APPENDIX A 


FINITE DIFFERENCE AND MATRIX ELEMENT EQUATIONS 
FOR THE NONSIMILAR BOUNDARY- LAYER EQUATIONS 


The following finite-difference equations were used to reduce the 
nonsimilar boundary- layer equations to a system of linear algebraic equations. 
The appropriate grid system from sketch (b) is reproduced here with redundant 
grid notation for convenience in writing the difference equations. 


n+1^4 


n — $ 
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f J 
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w < 
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d ( 
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w 1 
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m m + 1 

Sketch (b) 
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4 (Ay) 2 
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G b G e 
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(All) 


GM Y = G b fu? [CM c ' M a) + % - M d)]j 

Equations (A2) through (A4) , the Crank-Nicholson formulas, are used to approxi- 
mate the unknowns and derivatives of the unknowns. The central difference 
equations, (A5) and (A6) , are used to approximate derivatives of the proper- 
ties. Equations (A7) through (All) were chosen so that the unknowns appear 
linearly. 

These difference equations are substituted into equations (9) through 
(11) to obtain 

Momentum: 


A nF d + B n F e + C n F f = D x (A12) 

Energy: 

A 22 T d + B 22 T e + c 22 T f + A 2 iF d + B 2 iF e + C 21 Ff + A 23 C ld + C 23 C lf = D 2 (A13) 

Species : 

A 33 C l d + B 33 C l e + C33 C lf = °3 (A14) 

Expressions for the coefficients are presented below. There are 3(N - 2) 
equations for N nodal points and these are conveniently expressed by the 
matrix equation 


An w n-l + B n w n + Cn^n+1 -B n n = 2, 3, . . . , N - 1 


(A15) 


where 
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or 
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The matrix elements are given in the following equations, where for conve- 
nience, Z replaces the term a(l - y) . Also it is to be noted that follow- 
ing the first iteration quantities not subscripted with b are evaluated at 
(n, m + 1/2) by averaging quantities at n, m with those at (n, m + 1) . 

An = ^' Z + yj - C-Z4> y + a<|> - f - 2Sf s )l 

l 2 Ay 2 4 A y J 




B 2 2 


4>Z 2 
Pr Ay 2 


2SF b \ 

AS J 


C 22 

A 21 


^-A 2 2 + 

- 4>u e 2 Z 2 CF b ) y -| 
2Cpt 0 Ay J 


»Z 2 \ 
Pr Ay 2 / 



23 



C 23 = ( _A 23) 
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APPENDIX B 


FINITE DIFFERENCE AND MATRIX ELEMENT EQUATIONS 
FOR THE SIMILAR BOUNDARY- LAYER EQUATIONS 


The finite-difference equations used to reduce the similarity boundary- 
layer equations to a system of linear algebraic equations are given below. 

The grid system is shown in sketch (c) with redundant notation for conve- 
nience in writing the difference equations. The superscript (v) is understood 
to be the iteration number. . 



Sketch (c) 
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Substituting the above difference equations into equations (9) through (11) 
with the right-hand sides set to zero and approximating f(M)y by fj J v (M)y 

results in equations of the form given by equations (A12) to (A14) . The 
expressions for matrix elements are 
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